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1  Objectives 

Feedback  control  systems  for  aerospace  applications  must  maintain  precise  control  despite  uncer¬ 
tain  operating  conditions  and  unanticipated  circumstances  such  as  battle  damage.  These  systems 
must  be  designed  to  perform  robustly,  despite  uncertain  design  models  and  difficult  to  analyze 
nonlinear  effects.  They  must  also  be  capable  of  learning  and  adapting  when  accumulating  data 
indicates  that  previous  models  must  be  abandoned  and  that  existing  control  strategies  must  be 
changed.  This  three-year  research  program  centered  on  developing  and  improving  the  engineering 
techniques  for  robust  control  design,  with  particular  emphasis  on  the  need  for  data-driven  design 
methods  well-suited  to  situations  in  which  available  mathematical  models  are  poor  or  unreliable. 
The  research  results  are  expected  to  facilitate  the  design  of  control  systems  that  learn,  discover 
and  evolve  in  order  to  compensate  for  the  effects  of  battle  damage,  equipment  failures  and  other 
changing  circumstances.  Potential  applications  include  aircraft  stability  augmentation  systems, 
highly  maneuverable  aircraft  design,  missile  guidance  systems,  and  precision  pointing  and  tracking 
systems. 

2  Executive  Summary 

Thirty-nine  publications  supported  under  AFOSR  Grant  F49620-98- 1-0026  have  either  appeared, 
been  submitted  or  are  currently  pending  publication  [l]-[39].  The  subject  matter  of  these  AFOSR 
supported  publications  may  be  roughly  grouped  as  follows: 

•  Robust  LMI/BMI/IQC  Multivariable  Control  Theory 

-  Controller  Synthesis  [2,  15] 

-  Uncertain  Time-delays  [9,  10,  21,  25,  29,  33,  35,  39] 

-  Nonlinear  Systems  [4,  12,  13,  25,  30,  31,  34,  38] 

•  Unfalsified-Control,  Learning,  Adaptation  &  Controller  Identification  [5,  8,  11,  16,  20,  22,  23, 
26,  27,  32,  36,  37] 

•  Control  Design  Applications  [1,  11,  20,  21,  32] 

•  Editorial,  Review  &  Survey  Articles  [3,  7,  6,  14,  17,  19,  24,  28] 

Generally,  the  theoretical  developments  embodied  in  the  above  listed  recent  AFOSR-supported 
publications  have  been  accompanied  by  software  implementation  and  test  case  studies.  The  LMI 
/BMI/IQC  theory  plays  a  critical  role  in  extending  and  generalizing  the  Hoo  robust  controller 
design  theory  that  has  already  proven  its  value  in  aircraft  flight  control  applications  [40,  41,  42].  It 
allows  greater  flexibility  in  handling  structured  uncertainties,  controller  complexity  constraints  and 
gain-scheduling  requirements.  The  generalized  Popov  multiplier  robustness  analysis  and  synthesis 
techniques  [43,  44,  45]  developed  in  have  led  directly  to  improved  approaches  for  the  design  of  active 
vibration  damping  systems  for  flexible  space  structures  [1].  The  effective  and  rapid  transition  from 
theory  to  practice  has  been  facilitated  by  my  on-going  non-AFOSR-supported  involvement  with 
Dr.  R.  Y.  Chiang  in  upgrading  the  Matlab  Robust  Control  Toolbox,  a  robust  control  design 
software  product  published  by  The  MathWorks  and  in  widespread  use  by  government,  university 
and  aerospace  engineering  company  labs  [46]. 
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2,1  Research  Highlights 
Nonlinear  Robustness  Analysis 

A  major  breakthrough  has  been  the  development  of  optimally  precise  multiplier  methods  for  an¬ 
alyzing  robustness  of  systems  having  repeated  monotone  nonlinearities,  such  as  saturations  or 
relays  [34,  13],  We  have  derived  the  largest  class  of  bounded  linear  time-invariant  operators  that 
preserve  positivity  of  repeated  monotone  nonlinearities.  It  reduces  conservativeness  of  IQC  stabil¬ 
ity  analysis  of  systems  having  repeated  monotone  nonlinearities  to  the  minimum  possible.  Utility 
of  this  characterization  stems  from  the  fact  that  such  nonlinearities  are  typically  encountered  in 
stability  of  anti-reset  windup  schemes  among  other  applications.  We  have  also  derived  the  time 
varying  counterpart  of  this  class  of  operators. 

These  results  dovetail  with  earlier  results  giving  a  multiplier  interpretation  of  robustness  theory 
in  general,  including  its  newly  popular  LMI/BMI/IQC  formulations.  Our  results  enable  accurate 
reliable  analysis  of  stability  robustness  of  systems  like  missiles,  aircraft  and  robots  that  typically 
have  control  actuators  or  other  components  with  position  or  rate  limits  (i.e.saturation)  or  on-off 
(i.e.,  relay)  devices.  These  results  are  surprising  because  (1)  they  are  a  quantum  leap  forward 
in  a  once  much-studied  field  that  had  lain  dormant  for  most  of  the  past  30  years  and  (2)  they 
will  have  an  very  broad  practical  impact  because  will  give  significantly  less  conservative  sufficient 
conditions  for  stability  for  the  most  commonly  encountered  sorts  of  nonlinearities.  The  results  are 
the  best  possible  for  repeated  monotone  nonlinearities  in  the  sense  that  there  is  no  broader  class 
of  multiplier- type  (or  IQC)  robustness  conditions  that  can  do  better  without  further  information 
about  he  nonlinearities. 

Data-Driven  Unfalsified  Control 

Perhaps  the  most  significant  new  results  pertain  to  our  recent  development  of  data-driven  methods 
for  robust  and  adaptive  control  design  based  on  our  unfalsified  control  theory  [5,  8,  11,  16,  20,  22, 
23,  26,  27,  32,  36,  37]. Unfalsified  control  theory  facilitates  the  representation  of  adaptive  processes 
of  control  law  discovery  from  evolving  information  fiows  and  noisy  data.  In  this  paper,  the  theory 
of  unfalsified  adaptive  control  is  examined  from  the  behavioral  perspective  of  Willems.  An  ab¬ 
stract,  but  parsimonious,  min-max  optimization  problem  formulation  is  developed  that  describes 
and  unifies  direct  adaptive  control,  learning  theory  and  system  identification  problems  in  a  common 
behavioral  setting  based  on  the  concept  of  controller /model  unfalsification.  Thus,  adaptive  control 
is  seen  to  be  firmly  and  directly  linked  to,  and  to  conceptually  unified  with,  the  growing  body  of 
knowledge  on  behavioral  approaches  to  model  validation  and  unfalsified  system  identification.  The 
results  elucidate  and  underscore  the  fertile  conceptual  links  that  exist  between  adaptive  control 
theory  and  the  rich  theory  of  system  identification. 

This  unfalsified  control  theory  is  a  precise,  data-driven  approach  to  adaptive  controller  synthesis 
based  on  evolving  measurements  of  plant  response.  Unlike  traditional  control  design  methods  where 
controller  choices  generally  depend  heavily  of  prior  knowledge  of  plant  models  and  error-bounds, 
the  unfalsified  control  theory  gives  primary  emphasis  to  the  precise  analysis  of  the  implications 
of  evolving  experimental  measurement  data.  A  plant  model,  though  useful,  is  not  essential  and 
common  adaptive  control  pitfalls  associated  with  umealistic  assumptions  and  excessive  reliance  on 
prior  knowledge  are  circumvented. 
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2.2  Personnel 

Personnel  supported  by  AFOSR  Grant  F49620-98-1-0026  included  the  Michael  G.  Safonov  (PI), 
and  PhD  students  V.  Kulkarni,  M.  Jun,  S.  Bohacek,  A.  Grishencko.  The  research  supported  by 
the  grant  resulted  in  two  PhD  theses,  including  V.  Kulkarni  [38]  and  M.  Jun  [39].  Three  other 
fellowship  students  worked  closely  with  the  PI  on  the  research,  including  H.  Meng,  P.  Brugarolas 
and  R.  Mancera  whose  PhD  theses  expected  to  be  completed  in  2002. 

3  Major  Accomplishments 

3.1  All  Multipliers  for  Repeated  Monotone  Nonlinearities 

Introduction 

In  stability  analysis,  a  given  system  S  is  often  decomposed  into  two  interconnected  subsystems 
—  a  linear  time  invariant  subsystem  H  in  the  feedforward  path  and  an  otherwise  subsystem  A  in 
the  feedback  path.  More  often  than  not,  a  repeated  monotone  nonlinearity,  say  N,  is  encountered 
as  the  subsystem  A  (see,  e.g.,  [47],  [48],  [49],  [50],  [51]  and  references  therein).  A  key  step  in 
multiplier  based  stability  analysis  of  such  systems  is  to  characterize  a  class  of  multipliers,  i.e.  a 
class  of  convolution  operators,  such  that  every  element  M  of  it  preserves  positivity  of  N  in  the  sense 
that  N  >0  implies  M*N  >  0.  Stability  of  the  system  is  then  deduced  if  there  exists  at  least  one 
such  multiplier  M  such  that  MH  >  0  and  if,  in  addition,  H  has  a  finite  gain  (see  [52,  Theorem  2], 
[47],  [4]  and  references  therein  for  a  detailed  relevant  discussion).  Effectively,  positivity  preserving 
multipliers  give  an  integral  quadratic  constraint  (IQC)  characterization  of  N  (see  [53]  for  IQC’s  — 
theory  and  applications) .  The  larger  the  class  of  the  positivity  preserving  multipliers  the  better  it 

is,  for  the  sharper  is  its  IQC  characterization  and  the  lesser  is  the  conservativeness  in  the  stability 
analysis. 

The  best  available  class  of  positivity  preserving  multipliers  so  far  for  repeated  SISO  monotone 
nonlinearities  is  the  one  recently  derived  by  D’Amato  et  al,  [47].  Whether  it  is  the  best  possible 
as  well  has  been  unclear.  It  turns  out  that  they  have  stipulated  an  unnecessary  condition  on  the 
multipliers.  Identifying  and  relaxing  this  condition,  in  this  note  we  have  obtained  a  larger  —  indeed, 
the  largest  possible  —  class  of  positivity  preserving  multipliers  for  such  nonlinearities.  Specifically, 
we  have  characterized  the  largest  possible  classes  of  both  linear  time- invariant  as  well  as  linear  time 
varying  operators  that  preserve  positivity  of  such  nonlinearities.  Essentially,  our  results  generalize 
the  non-repeated  monotone  nonlinearity  results  of  Willems  [54,  Ch.  3]  to  the  case  of  repeated 
monotone  nonlinearities. 

Saturation  nonlinearities,  dead  zone  nonlineanities,  sigmoidal  nonlinearities  are  some  of  the 
many  examples  of  monotone  nonlinearities.  When  input-output  channels  of  A,  or  a  sub  block  of 

it,  feature  the  same  such  nonlinearity,  an  instance  of  repeated  monotone  nonlinearity  is  on  hand. 
Computation  of  stability  margin  of  anti  windup  schemes  is  one  of  the  engineering  applications  in 
which  repeated  monotone  nonlinearities  appear  (see,  e.g.,  [47,  49,  55]).  Reduction  of  conservatism 
in  such  stability  margin  estimates  is  thus  a  motivating  application  of  this  paper. 

The  paper  is  organized  as  follows.  In  Section  3.1,  the  necessary  terminology  is  introduced  and 
the  problems  are  formally  posed  in  Section  3.1.  Background  results  are  in  Section  3.1.  Our  main 
results  are  presented  in  Section  3.1  and  discussed  in  Section  3.1.  The  paper  is  concluded  in  Section 
3.1. 
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Preliminaries 

The  notation  used  is  summarized  in  Table  1.  Capital  letter  symbols,  e.g.  F  and  G,  denote  operators 
whereas  small  letters,  e.g.  x  and  y,  denote  real  signals  which  may  possibly  be  vector-valued  or 
matrix- valued.  The  vector  space  1.2  i®  generally  referred  to  as  £2-  Z  denotes  the  set  of  all  integers. 
(•)*  denotes  conjugate  transpose  of  a  vector  or  matrix  (•);  (■)^  denotes  its  transpose.  A  sequence 
{^(^)}^-oo  described  simply  as  {x}.  Statements  of  the  form  “A  related  to  B”  and  “C  related 
to  D”  are  abbreviated  as  “A  (C)  related  to  B  (D)” .  D  €  is  said  to  be  Hurwitz  if  each  of  its 
eigenvalues  has  a  strictly  negative  real  part.  Other  terms  not  defined  here  may  be  found  in  [56] 
and  [53]. 


Symbol 

ir(c) 

Z 

herm(m) 

skew(m) 


{^,y) 


£2 

T2 

x{-) 

S 

MD] 

MIMO 

SISO 


Meaning _ Table  1:  Notation 

Set  of  all  real  (complex)  numbers. 

Set  of  all  integers. 

=  Um  +  m*),  for  m  £  or 

=  |(m  -  m*),  for  m  €  or 

/  00 

y~]  y(k)^ x{k)  for  discrete  time  signals; 


fc  =  — 00 


/OO 

yitf  x{t) 

•OO 

=  \/{x,x). 

OO 


dt  for  continuous  time  signals. 


if  X  discrete  time; 


if  X  continuous  time. 


Space  of  discrete  time  signals  x  for  which  ||a;||  exists. 
Space  of  continuous  time  signals  x  for  which  ||a:||  exists. 
Fourier  transform  of  x,  either  discrete  or  continuous. 

_  J  Kronecker  S{k),  if  discrete  time; 

(  Dirac  5(t),  if  continuous  time, 

j-th  eigenvalue  of  matrix  H. 

Least  eigenvalue  of  a  Hermitian  matrix  D. 
Multi-Input-Multi-Output. 

Single-Input-Single-Output. 


Definition  1  [operator:  positive,  bounded] 

An  operator  F  mapping  a  space  X  into  itself  is  said  to  be  positive  if  {x,  Fx)  >  0  Vx  G  A.  A  set 
S  is  said  to  be  bounded  if  there  exists  a  7  G  such  that  ||y[|  <  7  for  all  y  E.  S.  An  operator 
F  :  X  -*Y  is  said  to  be  bounded  if  the  image  under  F  of  every  bounded  subset  of  X  is  a  bounded 
subset  ofY.  □ 

Definition  2  [sequences:  similarly  ordered,  unbiased] 

The  sequences  {a:}  and  {y}  of  real  scalars  are  said  to  be  similarly  ordered  if  x{k)  <  x{l)  implies 
y{k)  <  y{l)  for  all  k,l  E  Z.  They  are  said  to  be  unbiased  if  x{k)y{k)  >  0  Vfc.  They  are  said  to 
be  similarly  ordered  and  symmetric  if  they  are  unbiased  and,  in  addition,  the  sequences  {[a;]}  and 
{[j/]}  are  similarly  ordered.  □ 
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Definition  3  [associated  matrix,  kernel] 

Given  a  bounded  possibly  time  varying  linear  operator  M  :  ^2  ^2>  V  ~  given  as 

00 

y(^)  =  x{i)  yk  G  z, 

l=—oo 

where  rnk,i  €  yk,l;  the  associated  matrix  M  of  M  is  defined  as 


m_i  _i 

m-ifi 

m-1,1 

m_i,2 

mo,-i 

rhofi 

mo,i 

mo, 2 

mi  _i 

mifi 

mifi 

Wl,3 

m2,_i 

m2,Q 

^2,1 

m2, 2 

m2, 3 

The  symbol  mij,  i,j  G  Z  denotes  the  (i,j)-th  scalar  element  of  the  matrix  M;  for  example,  moo 
denotes  the  upper  left  entry  in  the  p  x  p  matrix  mo,o  o,nd  m_p^o  denotes  the  upper  left  entry  in  the 
p  X  p  matrix  m_i,o.  If  ink, l  =  ink+n,i+n  yk,l,n  G  Z  then  M  is  said  to  be  block  Toeplitz  and  M 
is  said  to  be  a  time  invariant  operator  or,  alternatively,  a  convolution  operator.  For  a  bounded 
possibly  time  varying  continuous  time  linear  operator  M  :  C2  C2 

/oo 

m{t,T)x(T)  dr  Vt  G  R. 

-00 

the  kernel  m{t,  r)  G  RP'^p  is  the  counterpart  ofrhk^i.  In  the  continuous  time  case,  M  is  called  a  time 
invariant  operator  or,  alternatively,  a  convolution  operator  ifm{t,  r)  =  m[t  +  v,t  +  v)  Vt,  t,u  sR. 
For  a  convolution  operator  M,  a  shorthand  notation  form{t,T)  and  rhij  is  rn{t  —  r)  and  Th{i  —j), 
respectively  with  rh{t)  and  rh{k)  denoting  the  respective  impulse  response.  □ 

Definition  4  [hyperdominance,  dominance] 

An  operator  M  ■.  I2  i2  is  said  to  be  doubly  dominant  if  the  elements  m^j  of  its  associated  matrix 
have  the  following  properties. 

00  00 

mu  >  53  -  53 

j=-oo,j#i  j=-oo,j^i 

If,  in  addition,  it  also  holds  that 

mij  <0,  yi^  j 

then  M  said  to  be  doubly  hyperdominant.  For  an  operator  M  :  C2  ^  C2,  these  notions  are  defined 
in  terms  of  its  kernel  in  an  analogous  manner  with  integrals  suitably  replacing  sums.  □ 


Definition  5  [monotone  nonlinearity] 

The  class  A/m  of  MIMO  monotone  nonlinearities  consists  of  all  memoryless  mappings  N  :  RP  -4  RF 
such  that  N{x)  is  the  gradient  of  some  convex  function  P  :  RP  — >  R  and  3C  G  RA  3  ||A/(x)||  < 
C||x||. 

M  =  {N  e  Mm\N{Q)  =  0},  Modd  =  {iV  G  N\N{x)  -  -iV(-x)  Vx}.  □ 
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Definition  6  [repeated  SISO  monotone  nonlinearity] 

The  class  of  repeated  SISO  monotone  nonlinearities  is  the  subclass  of  Af  with  element  N  G 

Af^^  of  the  form 

Nic)  =  mi)  HC2) . . .  mp)f  VC  G  (1) 

where  (p  G  Af,  4>  SISO.  A  shorthand  notation  for  (1)  is  N  —  diag{(l)).  The  class  Af^^  is  defined  by 
replacing  Af  in  the  definition  ofAf^^  by  Afodd-  D 


Definition  7  [multipliers] 

denotes  the  class  of  MIMO  convolution  operators,  either  continuous  or  discrete,  such  that 
the  impulse  response  of  an  M  £  AA^^  is  of  the  form 


m  =  g  5  —  h 

where  g,  h{-)  G  satisfy 

n  n 

(2) 

gu  ^  IffijI  + 

issl 

n  n 

II 

> 

(3) 

9n^  \9ij  \  +  II Ill 

The  subclass  AA^^  is  obtained  by  further  stipulating 

II 

(4) 

9ij  <  0  Vi  ^  j,  hij{-) 

Under  the  restriction 

IV 

0 

< 

(5) 

g,  h  are  Hermitian  matrices, 

(6) 

the  subclass  AA^  i-^odd)  derived  from  AA^  i-^odd)- 


Remark  1  D’Amato  et  al.  [47]  showed  that  AA^  (-^Sd)  preserves  positivity  of  Af^^  ■ 

Problem  Formulation 

Problem  1  Find  the  largest  class  of  bounded  linear  operators  and  the  largest  class  of  bounded 
convolution  operators  that  preserve  positivity  of  every  nonlinearity  in  Af^^ 

Background  Results 

Paraphrased  for  notational  ease,  the  main  result  of  D’Amato  et  ah,  viz.  [47,  Theorem  1],  is  as 
follows. 

Lemma  1  [47,  D’Amato  et  al.] 

AA^  (A4^d)  is  positivity  preserving  for  Af^^  (-^jdd)-  ^ 
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The  above  result  is  stated  a  sufficiency  condition  and  it  is  not  made  clear  if  there  exists  a  larger 
class  of  bounded  convolution  operators  that  preserves  positivity  of  this  regard,  it 

is  worthwhile  to  note  the  following  interesting  SISO  case  result  of  Willems  (see  [54,  Theorem  3.11, 
pp.  63]).  For  easy  reading,  its  statement  is  slightly  modified. 

Lemma  2  [54,  Willems] 

Let  M  :  £2  ^  I2  be  a  bounded  linear  operator.  Then,  {x,  My)  is  nonnegative  for  all  similarly  ordered 
unbiased  (similarly  ordered  symmetric  unbiased)  sequences  {m},  {y}  G  £2  if  and  only  if  M  is  doubly 
hyperdominant  (doubly  dominant).  □ 

Main  Result 

Theorem  1  [Solution  to  Problem  1] 

A  bounded  linear  operator  M  mapping  £2  into  £2  [or  £2  into  C2]  preserves  positivity  of  every 
N  G  {N  G  if  ii^  associated  matrix  [kernel]  is  doubly  hyperdominant  (doubly 

dominant).  Furthermore,  a  bounded  convolution  operator  M  mapping  £2  into  £2,  or  mapping  £2 
into  £2,  preserves  positivity  of  every  N  G  Af^^  {N  G  if  and  only  if  M  e  A4^^  (M  G  M^dd)- 

□ 


Proof:  We  shall  prove  the  result  for  Af^^.  The  case  for  Af^d  fohows  on  similar  lines.  First, 

the  result  will  be  proved  for  the  discrete  time  case. 

An  N  e  A/^^  can  be  expressed  as  IV  =  diag(</>)  where  ([  £  Af,  (f  SISO.  Given  sequences  {xj} 
of  real  valued  scalars,  define  yi  =  (f){xi)  i  =  l,2,...,p.  Note  that  the  sequences  {xj}  and  {yi}  are 
similarly  ordered  and  unbiased  for  all  i  since  N  e  Af.  Define  x{k)  =  [xi(fc)  X2{k)  ...  Xp(fc)]^, 
y{k)  =  [yi(fc)  y2{k)  . . .  ypik)]^  for  all  fc  G  Z.  Note  that  the  sequences  {x}  and  are  similarly 
ordered  and  unbiased.  Observe  that  {x,My)  =  (x,My^  where  the  sequences  {x}  and  {y}  are 
defined  by 

x(fc)  =  [xi{k)  X2{k)  ...  Xp{k)f,  y{k)  =  [yi(A;)  y2{k)  ...  y^ik)]^  Vfc 


and  M  :  £2  £2  with  its  associated  matrix  same  as  the  associated  matrix  M  oi  M.  Since  M 

is  bounded,  M  is  a  bounded  operator  as  well.  By  Lemma  2,  (t,  My'^  is  nonnegative  if  and  only 

if  M  is  doubly  hyperdominant.  This  is  proves  the  result  for  bounded  linear  operators.  To  prove 
the  result  for  bounded  convolution  operators,  note  that  the  associated  matrix  M  of  a  bounded 
convolution  operator  M  is  block  Toeplitz.  Since  M  is  block  Toeplitz,  the  double  hyperdominance 
conditions  need  only  be  checked  on  a  block  of  p  columns  and  on  a  block  of  p  rows.  The  conditions 
are  precisely  the  ones  given  by  (3)- (5). 

To  prove  the  result  in  the  continuous  time  case,  note  that  continuous  time  signals  Xc,  yc  can 
be  sampled  with  sampling  interval  e  to  produce  discrete  time  signals  in  £2,  say  Xd^e>yd,e  such 
that  Xd^e{i^)  =  y/e  Xc{ke),  yd,e{k)  —  \/e  yd^e).  Likewise,  given  a  continuous  time  linear  operator 


Me  :  £2  — ^  £2, 


mc{t,r)y{T)  dr  Vt, 


may  be  discretized  as  Zd^e  =  Md^cyd,e  i-C-  as 

CX> 

zdA^^)  =  X]  ^k,iyd,dl)  VA; 

l=—oo 
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Figure  1:  The  feedback  system  S.  H  is  stable,  causal  and  linear  time  invariant;  N  G 
where,  for  all  k  and  I, 

rriki  =  -  /  /  rricit,  r)  dr  dt. 

^  l)e,fce]  J T£((l—l)t,lc] 

With  this  discretization,  the  continuous  time  inner-product  {xc,  McUc)  is  then  recoverable  as  the 
limit 

{Xc,  McVc)  =  lim  {Xd,e,  Md^eVd^e)  ■ 

€—*0 

Taking  the  limit  as  e  — >  0,  the  continuous  time  case  proof  then  follows  using  the  discrete  time  case 
arguments.  QED. 


Remark  2  By  Theorem  1,  M^fd  preserves  positivity  of  the  identity  matrix.  Thus,  every  element 
M  of  has  the  property  that  herm(m(-))  >0.  ■ 

Remark  3  From  Theorem  1,  it  follows  that  every  convolution  operator  that  preserves  positivity 
of  and,  in  addition,  has  a  Hermitian  frequency  response  matrix  is  an  element  of 

AA^  However,  A4^  i-^^dd)  strictly  a  subset  of  {■^^dd)  because  the  condition  (6) 

stipulated  by  D ’Amato  et  al.  [47]  is  unnecessary  for  the  positivity  preservation.  ■ 

Discussion 

By  letting  go  of  the  unnecessary  condition  (6)  stipulated  by  D’ Amato  et  al.  [47],  we  have  obtained  a 
larger  —  indeed,  the  largest  possible  —  class  A4^^  i-^odd)  multipliers  that  preserve  positivity  of 
AfRS  However,  an  incremental  improvement  can  be  claimed  only  if  it  can  be  demonstrated 

that  use  of  this  larger  class  of  multipliers  leads  to  a  further  reduction  in  conservativeness  of  IQC 
stability  analysis.  The  following  example  demonstrates  that  the  conservatism  is  indeed  further 
reduced. 


Example  1  Consider  the  feedback  system  S  (see  Fig.  1)  in  which  N  G  AA^^  and  H  is  the  trivial 
stable  memoryless  linear  operator  having  constant  frequency  response 


h{u)  = 


0.5  -1 

2  -0.6 


Vw. 


The  objective  is  to  determine  if  the  system  is  stable.  Stability  is  established  if  (cf.  [52,  Theorem 
2],  [53,  Theorem  1])  there  exists  a  positivity  preserving  multiplier  M  such  that  the  operator  MH 
is  positive,  i.e. 


herm  >  0  Vw. 


(7) 
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It  may  be  checked  that  —h{uj)  G  is  not  Hurwitz  so  that,  by  Lyapunov  Criterion  (see,  e.g., 

[57,  Theorem  2.6-1]),  herm  (^Ph{uj)j  ':f>  0  for  any  P  =  >  0.  In  view  of  the  condition  (7),  it  then 

follows  that  stability  cannothe  determined  using  M  G  since  (see  Remark  2)  m{u>)  —  m{uj)^  >  0 
for  every  M  ^  .  On  the  other  hand,  choosing  the  asymmetric  multiplier 


M  = 


1  0 
-0.9  1 


(8) 


which  may  be  verified  to  be  in  (so  that  M  G  as  well),  the  condition  (7)  is  satisfied  so 

that  stability  of  the  system  is  established  using  the  class  derived  in  this  note.  Even  if  the 

nonlinearity  N  were  odd,  using  similar  arguments  it  follows  that  stability  cannot  be  determined 
using  M  G  whereas  it  can  be  determined  using  the  multiplier  given  by  (8).  This  demonstrates 
an  example  in  which  l®3,ds  to  a  strictly  less  conservative  stability  analysis  than 


Remark  4  The  above  example  shows  that  our  main  result  incrementally  improves  upon  the  main 
result  in  [47].  Its  utility  can  be  seen  via  the  same  engineering  application  considered  by  D’Amato 
et  al.  [47,  Automatica  version] .  ■ 


Summary:  All  Multipliers  for  Repeated  Nonlinearities 

We  have  derived  the  largest  possible  class  of  bounded  MIMO  linear  operators  that  preserve  posi¬ 
tivity  of  repeated  monotone  nonlinearities.  Also  derived  are  the  largest  possible  classes  and 
of  bounded  MIMO  convolution  operators  that  preserve  positivity  of  repeated  monotone  and 
repeated  odd  monotone  nonlinearities,  respectively.  It  follows  that  (see  Remark  2)  modulo  the 
restriction  that  multipliers  have  a  Hermitian  frequency  response  for  all  frequencies,  D’Amato  et  al. 
[47]  have  actually  derived  the  largest  possible  classes  of  multipliers  that  preserve  positivity  of  such 
nonlinearities.  The  less  restrictive  nature  of  our  multipliers  produces  less  conservative  stability 
results.  We  have  demonstrated  this  conservativeness  reduction  via  an  example. 


3.2  Unfalsified  Control  —  A  Behavioral  Approach  to  Learning  and  Adaptation 
Introduction 

A  fundamental  activity  in  processes  of  system  identification  and,  more  generally,  in  all  processes 
of  scientific  discovery  is  the  fitting  of  models  to  data.  Goals  or  cost  functions  are  set  based  on 
beliefs  about  instrument  accuracy  or  other  model  performance  criteria,  then  and  scientists  attempt 
experimental  validation,  or  unfalsification,  of  data  against  various  parameterized  classes  of  plausible 
models  in  the  hope  the  one  or  more  of  the  hypothesized  plausible  models  proves  have  a  superior  fit 
to  the  data,  relative  to  the  specified  modeling  goals  or  cost  functions.  The  challenges  faced  by  the 
to  the  system  identification  specialist  and  the  experimental  scientist  are  the  same.  Moreover,  these 
challenges  are  not  dissimilar  to  those  faced  by  adaptive  control  designers  seeking  to  find  a  faithful 
mathematical  representation  of  the  control-decision-relevant  information  in  evolving  observational 
data. 

One  interesting  development  in  recent  years  has  been  the  advent  of  the  unfalsified  control  para¬ 
digm  [58]- [5]  which  has  advanced  the  model  validation/unfalsification  paradigm  of  system  identifi¬ 
cation  theory  to  the  control  validation  paradigm  for  understanding  and  analyzing  adaptive  control 
algorithms.  In  adaptive  control  and  system  identification,  as  in  other  scientific  endeavors,  a  parsi¬ 
monious  mathematical  representation  of  the  essential  issues  is  preferred.  For  adaptive  control,  one 
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such  paradigm  is  provided  by  unfalsified  control  theory.  The  unfalsified  control  theory  views  the 
control  problem  as  an  identification  problem  in  which  the  objective  is  that  of  directly  identifying 
a  control  law  or  “action  rule”  that  is  consistent  with  traditional  control  performance  goals,  prior 
knowledge,  and  evolving  observational  data. 

Unfalsified  adaptive  control  software  has  been  developed  and  design  studies  have  been  con¬ 
ducted.  For  example,  the  theory  has  recently  been  applied  to  the  design  of  a  robust  adaptive 
controller  for  a  two-link  robot  manipulator  arm  in  [11]  and  to  the  design  of  an  adaptive  missile 
autopilot  [20],  Unfalsified  control  theory  provided  the  basis  for  a  general  purpose  algorithm  for 
automatic  tuning  of  PID  controller  gains  [22].  Recently,  the  unfalsified  control  approach  has  been 
applied  in  experimental  settings  ranging  from  the  ACC  benchmark  control  problem  [62]  to  indus¬ 
trial  process  control  [64,  65,  66].  These  design  studies  seem  to  confirm  theoretical  expectations  that 
adaptive  controllers  optimally  designed  via  unfalsified  control  theory  exhibit  a  precise,  sure-footed 
response  in  the  face  of  evolving  uncertainties  and  parameter  variations. 

In  this  paper,  we  formulate  the  unfalsified  adaptive  control  problem  in  the  behavioral  framework 
of  Willems  ([67], [68]).  Our  results  are  closely  related  to,  but  different  from,  the  indirect  adaptive 
method  of  Polderman  [69]  in  which  an  unfalsified  plant  model  is  identified  from  within  a  prescribed 
model  set.  Unlike  Polderman,  we  bypass  the  intermediate  step  of  plant  model  identification  and, 
also  unlike  Polderman,  we  make  no  assumptions  about  the  ‘true  plant’  lying  in  an  assumed  model 
set.  We  make  no  assumptions  on  the  plant. 

Background:  Behavioral  Theory 

At  the  heart  of  the  behavioral  theory  of  Willems  [67,  68]  is  the  definition  of  a  mathematical  model. 
This  definition  is  formulated  according  to  the  black  box  point  of  view,  “in  which  we  focus  on  how 
a  system  behaves,  on  the  way  it  interacts  with  its  environment,  instead  of  trying  to  understand, 
in  the  tradition  of  physics,  how  it  is  put  together  and  how  its  components  work”  ([67], [68])  .  This 
definition  of  a  mathematical  model  formalizes  the  black  box  point  of  view.  Like  Willems,  we  back 
off  “from  the  usual  input/output  setting,  from  the  processor  point  of  view,  in  which  systems  are 
seen  as  influenced  by  inputs,  acting  as  causes,  and  producing  outputs  through  these  inputs,  the 
internal  conditions,  and  the  system  dynamics.” 

Willems  begins  with  the  assumption  that  there  is  a  phenomenon  to  be  modeled.  He  then  “casts 
the  situation  in  the  language  of  mathematics  by  assuming  that  the  phenomenon  produces  elements 
in  a  set  Z”  ([67], [68]),  called  the  universum.  The  elements  of  Z  are  called  the  outcomes  of  the 
phenomenon.  “A  (deterministic)  mathematical  model  for  the  phenomenon  (viewed  purely  from  the 
behavioral,  the  black  box  point  of  view)  claims  that  certain  outcomes  are  possible,  while  others  are 
not.  Hence  a  model  recognizes  a  certain  subset  B  of  Z.  This  subset  will  be  called  the  behavior  (of 
the  model).”  Formally, 

Definition  3.1  A  mathematical  model  is  a  pair  (Z,B),  with  Z  the  universum  —  its  elements  are 
called  outcomes  —  and  B  CZ  the  behavior. 

Definition  3.2  A  controller  is  a  mathematical  model. 

Regarding  data  and  measurements,  Willems  [67]  says;  “We  will  now  cast  measurements  in  this 
setting.  We  will  assume  that  we  make  certain  measurements  which  we  will  call  the  data.”  “. .  .we 
. . .  assume  that  the  data  consists  of  observed  realizations  of  the  phenomenon  itself.  Thus,  a  data 
set  will  be  a  nonempty  subset  X>  of  Z.  ” 
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Following  Safonov  and  Tsao  [58] ,  we  will  work  with  data  information  that  can  evolve  with  time. 
Thus  we  will  have  a  universum  of  time  signals  and  a  data  set  contained  in  a  time  varying  projection 
of  this  universum. 

Definition  3.3  Given  a  vector  space  of  time  signals  Z,  a  model  a  mapping  Pt  :  Z  Z 

and  a  data  set  Vr  =  PrP  C  Pt(Z.),  we  say  that  the  model  {Z,B)  is  unfalsified  by  the  data  set  Vr  if 

Vr  C  Pr{B). 

Typically  Pt{x)  is  the  experimental  observation  time  sampling  operator,  which  returns  values  of 
x{t)  only  for  past  time  instants  (or  possibly  time  intervals)  over  which  experimental  observations  of 
x{t)  have  been  recorded.  In  this  setting,  definition  of  controller  falsification  (cf.  [58,  60])  becomes 

Definition  3.4  Given  a  vector  space  of  time  signals  Z,  a  controller  (Z,Bc),  a  desired  closed  loop 
behavior  (Z,Pd),  a  mapping  Pr  :  Z  Z,  and  a  data  set  Vr  C  Pt{Z.),  we  say  that  a  controller 
(Z,  Be)  is  unfalsified  by  the  data  set  Vr  if 

Pri{PrHT^r))r)Be)cPr{Bd). 

The  data  set  Vr  is  a  set  of  actual  experimental  observations  of  the  plant  behavior  as  observed 
through  the  time-sampler  Pr-  Thus,  Pf^{Vr)  is  the  set  of  behaviors  that  interpolate  the  observed 
data.  For  example,  if  we  have  recorded  experimental  observations  of  the  first  component  xi{t)  of  a 
vector-valued  signal  x{t)  =  [xi{t) , X2{t) , . . . , Xn{t)]'^  €  L2[0,oo)  during  the  time  interval  t  G  [0,5], 
then  Pf^{x)  is  the  set  of  signals  {y  G  L2[0,oo)  j  yi{t)  =  xi{t)Wt  G  [0,5]}.  The  set  Be  is  the  set  of 
signals  which  satisfy  the  constraints  imposed  by  the  controller  c,  so  definition  3.4  says  roughly  that 
a  controller  is  defined  to  be  unfalsified  if  the  set  of  signals  x  that  are  consistent  with  the  data  and 
the  controller  is,  at  the  past  observation  times,  a  subset  of  a  given  performance  target  set  Pr{Bd). 

A  particularly  useful  projection  operator  for  dealing  with  past  time  only  information  is  the  time 
truncation  operator  Pr  defined  by 

As  explained  by  Willems  ([67], [68]),  the  intersection  of  behaviors  is  “a  way  of  formalizing  that  addi¬ 
tional  laws  are  imposed  on  a  system.”  Thus,  the  role  of  a  controller  is  to  impose  constraints  on  the 
plant  behavior.  On  the  other  hand,  our  goal  is  to  select,  based  on  the  data,  the  constraints  imposed 
by  the  control  law  and  the  performance  criterion,  the  best  among  the  set  of  given  controllers.  In 
order  to  do  that  we  introduce  a  cost  function 

J{z)  :  Z  (9) 

which  may  be  used  to  sift  controllers  and  to  choose  an  optimal  cost-minimizing  controller  having 
the  least  unfalsified  cost  based  on  the  experimental  evidence  Vr- 

Direct  Adaptive  Control:  Behavioral  Formulation 

We  now  explain  how  problems  of  adaptive  control  and  learning  theory  may  be  parsimoniously  and 
precisely  embedded  within  the  behavioral  framework.  At  the  crux  is  the  observation  that  most 
such  problems  may  be  faithfully  represented  in  terms  of  constraints  on  signals  and  other  variables 
z  £  Z.  The  set  Z  is  called  the  universum.  Typically,  the  n-tuple  z  includes  directly  observable 
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manifest  variables  ^manifest  (viz.,  control  and  sensor  signals),  command  input  signals  ^command,  and 
possibly  additional  latent  variables  Z|ate„t  such  as  disturbances,  noise,  state-variable  trajectories, 
error  signals  and  so  forth.  That  is,  Z  —  { '^manifest i '^command)  ^latent}  ^  2[n,anifest  X  Z,,Qmniand  ^  ^latent  1^ 

stochastic  settings,  the  latent  variable  n-tuple  Zi^tent  also  includes  conditional  probability  density 
functions  describing  some  of  the  other  latent  variables  [60,  18,  8] . 

It  is  convenient  to  view  the  constraints  on  the  n-tuple  of  signals  2  G  Z  as  arising  from  four 
distinct  types  of  information,  each  possibly  evolving  with  time: 

1.  Goal  (cost  function  and/or  design  specification) 

2.  Belief  (assumptions,  prior  knowledge,  noise  models,  plant  parameterizations,  etc.) 

3.  Hypothesis  (candidate  control  law,  candidate  plant/noise  model) 

4.  Data  (observations,  samples  of  the  signal  available  at  current  time  r) 

Each  of  these  four  types  of  information  is  representable  as  a  mathematical  constraint  on  z: 
V^command,  J (z)  <7,  (cost  J {z)  >  0  uo  bigger  than  7  for  any  command  input)  (10) 

K{z)  =  Q,  (hypothetical  controller  and/or  model  X  G  K)  (11) 

B{z)  <0,  (fixed  beliefs,  assumptions  &  prior  knowledge)  (12) 

-Pr(2manifest)  =  -^data,  (-^rranifet  Hiust  interpolate  observed  data  z^^ta)  •  (13) 

In  turn,  the  four  constraints  (10)-(13)  define,  respectively,  four  subsets  of  the  universum  Z,  viz. 

Zgoal(7))  ^  hypothesis  ^belief?  ^data  ^  Z.  (14) 

Thus,  the  problem  of  direct  adaptive  control  (or,  controller  identification  from  data),  can  be  for¬ 
mulated  in  the  Willems  behavioral  framework  as  follows. 


Problem  1  (Behavioral  Adaptive  Control)  Given  a  class  of  controllers  K  =  {{Z,Bc{d))  \  0  G 
©},  where  0  is  a  set  of  parameter  vectors,  the  performance  (cost)  index  J{z)  the  time  truncation 
operator  P-r,  r,  and  a  data  set  T>r  C  PrZ,  find  the  set  of  parameters  ©*  such  that  A  G  K  that 
minimizes  the  cost  7  subject  to  the  constraint  (cf.  [60,  18,  8])  that,  for  each 

C  ^  21hypothesis(A')  r*)  Z|,e|jef  ri  Z(j3,3,  (15) 

there  is  at  least  one  z  such  that 


and 


‘2'command  —  Scommand 

z  G  Zg(,3i(7)  n  2hypothesis(A'(0))  G  Z^eij^f  [7  Z^3,3. 


(16) 

(17) 


Discussion 

If  the  set  (15)  is  empty  for  some  K,  then  the  currently  available  data  Zdata  provides  no  information 
on  this  K,  which  is  therefore  trivially  optimal  with  cost  7  =  0;  otherwise,  the  adaptive  feedback 
control  problem  emerges  as  the  following  optimization:  At  each  time  t,  find  a  control  law  K  which 
solves 

7™,  :=  min  max  min  7  (18) 

K  ^  z 

subject  to  (15)-(17).  In  many  practical  cases  (e.g.,  [58,  18,  20,  11]),  the  cost  7  can  be  expressed 
directly  in  terms  of  Zd3,3  and  K  in  which  case  (15)-(18)  simplify  to  7opt  := 

Noteworthy  are  the  symmetries  revealed  in  the  condition  (17)  with  respect  the  information 
content  of  goal,  belief,  hypothesis  and  data.  Set  intersection  is  a  commutative  and  associative 
operation;  so  all  four  types  of  information  are  logically  equivalent  in  (17).  For  example,  this  means 
that  the  prejudice  inherent  in  viewing  one’s  data  through  a  prism  of  belief  is  logically  equiv¬ 
alent  to  assuming  additional  data  “interpolation”  constraints  (Zd3,3  <—  Z^^  ©  Zteiief).  The  standard 
unfalsified  control  problem  considered  in  [58,  59,  60,  20,  11,  5]  corresponds  to  the  limiting  case  in 
which  “the  prism  of  belief”  Zteiicf  is  the  unconstraining  “all-pass”  filter  Z  (i.e.,  the  universum). 


4  CONCLUSIONS 


16 


Equations  (15)-(17)  underscore  fertile  conceptual  links  between  adaptive  control  theory  and 
the  rich  theory  of  system  identification:  The  chief  differences  between  identification  and  adaptive 
control  arise  from  the  precise  forms  of  the  cost  functions  J{z)  and  of  the  admissible  hypotheses 
K{z).  In  system  identification  the  admissible  K{zys  are  typically  noisy  open-loop  plant  models 
and  the  cost  function  J{z)  measures  probable  modeling  error  deduced  via  a  prism  of  beliefs  about 
noise  statistics.  In  adaptive  control  on  the  other  hand,  the  admissible  ii'(z)’s  might  typically  be 
candidate  controllers  and  the  cost  J {z)  could  be  a  weighted  sum  of  the  sizes  of  tracking  error  and 
control  signals. 

Summary:  Unfalsified  Control  —  A  Behavioral  Approach  to  Learning  and  Adaptation 

The  main  goal  of  unfalsified  control  theory  has  been  to  close  the  loop  on  the  adaptive  and  robust 
control  design  processes  by  developing  data-driven  methods  to  complement  traditional  model-based 
methods  for  the  design  of  robust  control  systems.  The  crux  of  the  unfalsified  control  theory  is  the 
observation  that  adaptive  control  is  from  a  behavioral  theory  perspective  essentially  equivalent  to 
system  identification.  In  this  paper,  we  have  developed  a  behavioral  formulation  of  the  problem  of 
direct  adaptive  control,  viz.,  the  problem  of  identifying  an  optimal  controller  that  is  unfalsified  by 
data  available  at  each  time  r  with  respect  to  the  least  value  of  a  cost  function.  The  main  result  is 
the  formulation  of  direct  adaptive  control  problems  provided  by  Problem  1 .  This  result  establishes 
a  firm  theoretical  link  between  Willems’  behavioral  framework  and  direct  adaptive  control  theory, 
expanding  known  links  to  model  validation,  unfalsified  system  identification  theory,  and  behavioral 
indirect  adaptive  control  approach  of  Polderman  [69]. 

4  Conclusions 

With  support  from  AFOSR  Grant  F49620-98-1-0026,  significant  progress  has  been  made  in  theory 
for  reliable  computation  of  robust  controllers,  in  the  field  of  nonlinear  robustness  analysis,  and  in  the 
development  of  the  unfalsified  control  theory  formulation  of  adaptation  and  learning  problems.  Our 
new  results  in  nonlinear  robustness  analysis  for  systems  with  repeated  nonlinearities  are  important 
breakthroughs  of  broad  applicability  to  most  all  control  designs.  They  will  enable  one  to  more 
accurately  predict  and  correct  stability  problems  that  commonly  occur  with  control  actuators  hit 
position  or  rate  limits. 

Progress  in  unfalsified  control  theory,  a  new  formulation  of  adaptive  control  problems  developed 
with  AFOSR  support,  gives  sharp  mathematical  representation  of  the  role  of  experimental  data  in 
identifying  robust  control  laws  and  provides  a  practical  technique  for  identifying  robust  controllers 
in  real-time  with  little  or  no  apriori  information.  This  theory  cuts  to  the  heart  of  a  long-standing 
problem  faced  by  control  engineers,  which  has  been  the  need  for  a  data-driven  theory  that  provides 
a  unified  basis  for  precisely  representing  and  exploiting  evolving  information  flows  from  models, 
noisy  data,  and  more.  The  theory  has  unveiled  important  links  between  robust  control,  adaptive 
control  and  artificial  intelligence.  It  is  a  conceptual  breakthrough  because  it  distills  the  mathe¬ 
matical  essence  of  control-oriented  learning  by  focusing  sharply  on  what  is,  and  is  not,  knowable 
from  experimental  data  and  by  challenging  both  the  need  and  the  appropriateness  of  a  number  of 
common  assumptions.  The  results  of  our  unfalsified  control  research  have  not  only  solidified  the 
theoretical  foundation  of  adaptive  control.  They  have  also  permitted  efficient  design  of  feedback 
control  systems  endowed  with  an  unprecedented  ability  to  accurately  assess  and  exploit  evolving 
real-time  information  flows  as  they  unfold,  thereby  endowing  control  systems  with  the  intelligence 
to  quickly  and  sure-footedly  adapt  to  unfamiliar  environments.  The  effectiveness  of  the  unfalsified 
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approach  to  adaptive  control  design  has  recently  been  demonstrated  with  simulated  design  studies 
for  missile  autopilot  [20],  a  robot  manipulator  arm  [11],  an  industrial  machine  tool  [66]  and  a  gen¬ 
eral  purpose  adaptive  PID  controller  that  automatically  discovers  suitable  controller  gains,  when 
they  exist,  based  on  real-time  data  without  the  need  for  plant  model  [22].  These  research  results 
will  mean  control  designs  that  can  more  quickly  and  reliably  compensate  for  uncertain  effects  of 
battle  damage,  equipment  failures  and  other  changing  circumstances. 
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Interactions  &  Coupling  Activities 

Gov’t  Panel:  Michael  G.  Safonov 

NSF  Panel  on  Knowledge  and  Distributed  Intelligence. 

Learning  and  Intelligent  Systems. 

Arlington,  VA,  July  1998. 

Talk:  M.  G.  Safonov. 

Robust  control,  feedback  and  learning. 

Talk,  Workshop  for  Michael  Athans,  Tampa,  FL,  December  19,  1998,  1998. 

60th  birthday  celebration  honoring  M.  Athans. 

Talk:  M.  G.  Safonov. 

Synthesis  of  positive  real  feedback  systems:  A  simple  derivation. 

Talk,  UCSB,  Santa,  Barbara,  CA,  October  16,  1998. 

Gov’t  Panel:  NSF  Panel.  Knowledge  and  Distributed  Intelligence:  Learning  and  Intelligent  Sys¬ 
tems,  Washington,  DC,  March  15-16,  1999. 

Talk:  M.  G.  Safonov. 

Robust  control,  feedback  and  learning. 

Talk,  IFAC  Inti.  Workshop  on  Control  of  Uncertain  Systems,  Hong  Kong,  Univ.  of  Science  k 
Technology,  June  30  -  July  2,  1999. 

Talk:  M.  G.  Safonov. 

Robust  control,  feedback  and  learning. 

Talk,  Poster  Presentation,  AFOSR  Workshop  on  Dynamics  and  Control,  Dayton,  OH,  August 
4-6,  1999. 

ftp : //routh . use . edu/pub/saf  onov/saf  o99h . ppt 

Talk:  M.  G.  Safonov. 

CACSD  design  process. 

Talk,  Panel  Discussion  on  Perspectives  on  Computer  Aided  Control  Systems  Design,  IEEE 
Symp.  on  Computer  Aided  Control  System  Design,  Kohala  Coast-Island  of  Hawaii,  HI,  Au¬ 
gust  22,  1999. 

Plenary  session. 
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Talk:  M.  G.  Safonov. 

Multiplier  IQC’s  for  uncertain  time  delays. 

Talk,  Tokyo  University,  Tokyo,  Japan,  May  22,  2000. 

Talk:  M.  G.  Safonov. 

Robust  control,  feedback  and  learning. 

Talk,  SICE  Conference,  Kariya,  Japan,  May  24,  2000. 

Invited  plenary  talk. 

Talk:  M.  G.  Safonov. 

Robust  control  tutorial. 

Talk,  Tokyo  University,  Tokyo,  Japan,  May  31,  2000. 

Talk:  M.  G.  Safonov. 

Zames-Falb  multipliers  for  MIMO  nonlinearities. 

Talk,  Kyoto  University,  Kyoto,  Japan,  June  2,  2000. 

Talk:  M.  G.  Safonov. 

Zames-Falb  multipliers  for  MIMO  nonlinearities. 

Talk,  Osaka  University,  Osaka,  Japan,  June  6,  2000. 

Talk:  M.  G.  Safonov. 

Unfalsified  direct  adaptive  control  of  a  two-link  robot  arm. 

Talk,  Titech,  Tokyo,  Japan,  June  13,  2000. 

Talk:  M.  G.  Safonov. 

Zames-Falb  multipliers  for  MIMO  nonlinearities. 

Talk,  AFOSR  Workshop  on  Dynamics  and  Control,  Dayton,  OH,  August  21-23,  2000. 

Talk:  M.  G.  Safonov. 

Robust  control,  feedback  and  learning. 

Talk,  University  of  California,  Santa  Barbara,  CA,  October  13,  2000,  2000. 

Talk:  M.  G.  Safonov. 

Robust  control,  feedback  and  learning. 

Talk,  Caltech,  Pasadena,  CA,  October  16,  2000,  2000. 

Talk:  M.  G.  Safonov. 

Zames-Falb  multipliers  for  MIMO  nonlinearities. 

Talk,  University  of  California,  Santa  Barbara,  CA,  November  17,  2000. 

Talk:  M.  G.  Safonov. 

Robust  control,  feedback  and  learning. 

Talk,  IEEE  Control  Society,  San  Diego  Chapter,  La  Jolla,  CA,  January  18,  2001. 

Talk:  M.  G.  Safonov. 

Zames-Falb  multipliers  for  MIMO  nonlinearities. 

Talk,  University  of  California  San  Diego,  La  Jolla,  CA,  January  18,  2001. 

Talk:  M.  G.  Safonov. 

Zames-Falb  multipliers  for  MIMO  nonlinearities. 

Talk,  Workshop  in  honor  of  Boyd  Pearson,  Rice  University,  Houston,  TX,  March  9-10,  2001. 
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Talk:  M.  G.  Safonov. 

Data-driven  behavioral  formulation  of  the  adaptive  feedback  control  problem. 

Talk,  AFOSR  Workshop  on  Dynamics  and  Control,  Dayton,  OH,  July  30  -August  02,  2001. 
ftp : //routh . use . edu/pub/saf onov/saf oOlo . pps 

Soft-ware:  Effective  and  rapid  transition  from  theory  to  practice  has  been  facilitated  by  my  on¬ 
going  non-AFOSR-supported  involvement  with  Dr.  R.  Y.  Chiang  in  upgrading  the  Matlab 
Robust  Control  Toolbox,  a  robust  control  design  software  product  published  by  The 
MathWorks  and  in  widespread  use  by  government,  university  and  aerospace  engineering  com¬ 
pany  labs. 

Technology  Transitions 

•  Semiconductor  Manufacturing  Process  Control 

Enabling  Unfalsified  Control  Theory 

Research: 

Performer:  Dr.  Robert  Kosut,  SC  Solutions,  Sunnyvale,  CA  94085;  (408) 

617-4527 

Customer:  Dr.  Robert  Kosut,  SC  Solutions,  Sunnyvale,  CA  94085;  (408) 

617-4527 

Result:  Used  unfalsified  control  technique  to  tune  controller  for  semi¬ 

conductor  manufacturing  process. 

Application:  Higher  quality  semiconductors. 

•  Commercial  Robust  Control  Software 

Enabling  Robust  control  theory,  including  LMI/BMI  methods 

Research: 

Performer:  Michael  G.  Safonov/USC  (213)  740-4455 

Customer:  The  MathWorks,  Inc./Natick,  MA  Mr.  John  Little,  (508) 

653-1415 

Result:  Improved  and  reorganized  algorithms  for  robust  control  in 

various  MATLAB  software  products,  including  the  Robust 
Control  Toolbox,  Mu-Synthesis  Toolbox,  and  LMI  Toolbox. 

Application:  Faster,  cheaper,  more  reliable  control  design  for  advanced 
aerospace  and  other  control  systems. 

New  Discoveries,  Inventions,  or  Patent  Disclosures.  None 

Honors  /  Awards 

•  Honor/ Award:  Elected  IEEE  Fellow  1989 

Honor/Award  Recipient(s):  Michael  G.  Safonov 

Awarding  Organization:  IEEE 

•  Honor/Award:  Author  of  4  papers  selected  for  inclusion  is  IEEE  1986 

Press  anthology  of  key  papers  in  control,  “Robust 
Control”  ed.  P.  Dorato 

Honor /Award  Recipient  (s):  Michael  G.  Safonov 
Awarding  Organization:  IEEE 
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•  Honor/Award:  Author  of  4  papers  selected  for  inclusion  is  IEEE  1990 

Press  anthology  of  key  papers  in  control,  “Recent  Ad¬ 
vances  in  Robust  Control”  ed.  P.  Dorato  and  R.  Ye- 
davalli 

Honor /Award  Recipient  (s):  Michael  G.  Safonov 
Awarding  Organization:  IEEE 

•  Honor/Award:  Awards  Chair  1993-1995 

Honor/ Award  Recipient  (s):  Michael  G.  Safonov 

Awarding  Organization:  American  Automatic  Control  Council 

•  Honor/Award:  Editorial  Board.  1990-2001 

Int.  J.  Robust  and  Nonlinear  Control 
Honor/ Award  Recipient  (s):  Michael  G.  Safonov 
Awarding  Organization:  J.  Wiley 

•  Honor/Award:  Assoc.  Editor,  Systems  and  Control  Letters  1995-^2001 

Honor /Award  Recipient  (s):  Michael  G.  Safonov 

Awarding  Organization:  Elsevier  Science  Publishers 


•  Honor/Award:  Invited  Plenary  Speaker  May  2000 

Honor/Award  Recipient(s):  Michael  G.  Safonov 

Awarding  Organization:  Society  of  Instrumentation  and  Control  Engineers  (SICE,  Japan). 


